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To unravel how the microstructure affects the fracture surface roughness in heterogeneous brittle 
solids like rocks or ceramics, we characterized the roughness statistics of post-mortem fracture 
surfaces in home-made materials of adjustable microstructure length-scale and porosity, obtained by 
sintering monodisperse polystyrene beads. Beyond the characteristic size of disorder, the roughness 
profiles are found to exhibit self-affine scaling features evolving with porosity. Starting from a null 
value and increasing the porosity, we quantitatively modify the self-affine scaling properties from 
anisotropic (at low porosity) to isotropic (for porosity larger than 10%). 

PACS numbers: 46.50.+a, 62.20.M-, 78.55.Qr 


I. INTRODUCTION 

Fractography, i.e. the morphological characterization 
of post-mortem fracture surfaces, is a classical tool used 
to identify the mechanisms and the damaging processes 
(fatigue, stress corrosion, cleavage, plastic cavitation, 
crazing,...) responsible of failure (see [1 1 for recent text¬ 
book). Since the pioneer work of Mandelbrot et al. j.2], 
numerous studies have evidenced the existence of self- 
affine scaling invariance over a wide range of length- 
scales (see mm for reviews). The early measurements in 
various materials including metallic alloys, ceramics and 
rocks EHZj reported values for the roughness exponent ( 
close to 0.8, suggesting the existence of a universal value 
13 El, independent of the propagation mode and material. 
These experimental observations yielded intense theoret¬ 
ical mm and numerical [T5HT8] researches. 

Recently, the picture has been made more complex. 
Firstly, the roughness exponent has been found to 
(slightly) depend on the anisotropy of the material mi¬ 
crostructure m and on the fracture speed [20]. Sec¬ 
ondly, fracture surfaces were shown to exhibit anisotropic 
scaling features, characterized by two different rough¬ 
ness exponents whether observed along the direction 
of crack front or along the direction of crack growth 
mm- Thirdly, fracture surfaces exhibit anomalous 
scaling [231 HU; the introduction of an additional global 
roughness exponent ( g i^ is then necessary to describe 
the scaling between the global crack width and the spec¬ 
imen size. Fourthly, multiaffinity, disappearing at large 
scale, was invoked m m- Finally, the scale invari¬ 
ance properties were found to depend on the propaga- 
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tion mode: The seminal self-affine feature with £ ~ 0.8 
is to be linked with quasi-brittle fracture; lower values 
for around 0.4 —0.5 [27H30| or even logarithmic rough¬ 
ness m are to be associated with brittle fracture; and 
a multi-affine regime with a roughness exponent close to 
0.5 is characteristic of ductile fracture [32U34] . Note that 
different regimes can be observed on a same fracture sur¬ 
face QUEUE!], depending on the scale of observation 
and its position with respect to the relevant sizes asso¬ 
ciated with the various dissipation mechanisms (plastic 
zone size, damage zone size, fracture process zone size). 
This enables to infer the fracture process zone from a 
statistical characterization of crack roughness mi- 


Most past studies dedicated to the roughness of cracks 
aimed at characterizing the scale-invariant properties of 
fracture surfaces. This quest of universality classes, 
i.e. of features independent of the fracturing conditions 
(loading and material parameters), strayed from the ini¬ 
tial metrology purpose of the fractography science. Here, 
we go back to this primary purpose and seek to charac¬ 
terize how the microstructure affects the fracture sur¬ 
face roughness. Rather than using existing materials like 
ceramics or rocks, we use home-made porous solids ob¬ 
tained by sintering spherical monodisperse grains. The 
advantage is that both the grain size and the porosity can 
be easily adjusted by modifying, respectively, the bead 
diameter and the sintering parameters, while presenting 
a structure of cemented grains very comparable to the 
structures found in natural heterogeneous brittle solids 
like rocks (sandstones for instance) or concrete, and other 
artificial heterogeneous brittle solids like pharmaceutical 
pills or sintered ceramics. The procedure for preparing 
the samples is described in detail in Sec |II A| The samples 
are finally broken using a technique detailled in Sec |IIB| 
that permits to grow stable Mo de I c racks. The frac¬ 
tography of the surfaces (See Sec IIC) is used to deter¬ 
mine the propagation mode of the fracture i.e. inter- 
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or intra-granular propagation. The spatial correlation of 
the roughness of the crack surfaces is, finally, character¬ 
ized by the structure function of profiles taken along and 
perpendicular to the direction of crack propagation. The 
evolution of this function as function of the bead size, the 
direction of measurements, and the porosity is reported 
in Sec |III| Finally, we discuss in Sec |IV| the change in 
the self-affine exponents observed when the porosity is 
modified. 


II. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

The experiments reported here consist (i) in driving a 
stable tensile crack with a controlled velocity in a mate¬ 
rial obtained by sintering polystyrene beads and (ii) sub¬ 
sequently, in analyzing the spatial distribution of crack 
roughness. The material fabrication, fracture set-up, 
and topographical recording of the postmortem rough¬ 
ness profiles are presented thereafter. 


A. Sintering protocol 

The samples are obtained by sintering monodisperse 
polystyrene beads (Dynoseeds©) of various diameters. 
Four steps are achieved: (i) an home-made mold whose 
geometry is given in Fig. & is filled with beads and 
heated up to 90% of the temperature at glass transition 
to soften the beads; (ii) a slowly linearly increasing com¬ 
pressive force is then applied up to the prescribed sinter¬ 
ing load F s for nearly one hour to achieve sintering; (iii) 
then, an annealing stage is performed: the mold is un¬ 
loaded, unscrewed and loosened, keeping the temperature 
high enough to avoid thermal-shocks; (iv) the sample is 
finally slowly cooled down to ambient temperature. The 
whole sintering protocol is sketched with more details in 
Fig. [lb, and the relevant temperatures are gathered in 
Tablejl]: (see also [38! for an extensive presentation). 

The bead diameter do was varied between 21 and 
583 /im. By changing F s between 0.1 and 8 kN (cor¬ 
responding to loading pressure ranging from 6 kPa to 
490 kPa), we are able to fabricate specimens of porosi¬ 
ties between ~ 0% and 19% (see Tab. [T|. Smaller values 
of F s leads to unusable friable samples. The porosity <f> is 
measured as the ratio = 1 — p/po of the density of the 
sample p and a reference density po corresponding to the 
density of the void free samples obtained with the highest 
sintering force. It is found that po is slightly bead size 
dependent, which may be explained by a slightly different 
initial bead material. 


B. Fracture wedge-splitting tests 

Specimens (Fig. Jlk) are parallelepipeds of dimension 
140 x 125 x IF mrmin the x xy x z direction, where the 
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FIG. 1. a) Geometry and dimensions of the wedge splitting 
samples, b) (Color online) Evolution sketch of the tempera¬ 
ture of the oven, of the sample, of the force applied during 
the sintering protocol (See [38] for details), c) Temperature 
value at the different stages of the sintering protocol for the 
different bead sizes. 


specimen width IF depends the material porosity and 
ranges from 15 mm (for the most compacted specimens) 
to 20 mm (for the most porous specimens). A rectangular 
notch is milled by cutting a 42 x 30 x IF mm 3 rectangle 
from the sample at one of the two lateral faces. A 10-mm- 
long 2-mm-thick groove is then introduced in the middle 
of the rectangular notch using a diamond saw and a seed 
pre-crack is added at the groove tip by means of a razor 
blade. This latter operation prevents the propagation of 
a dynamic fracture and allows to grow slow stable cracks. 
Thereafter, we define the axes x parallel to the direction 
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Experiment n° 

1 

2 

3 

4 

5 

6 

d 0 (/mi) 

21 

21 

42 

42 

42 

81 

F s (kN) 

4 

8 

1 

8 

8 

2 

V wdg (nm.s -1 ) 

16 

16 

16 

16 

16 

16 

$ (%) 

1 

0 

1 

0 

0 

1 

Experiment n° 

r 

8 

9 

10 

11 

12 

d 0 (/mi)l 

81 

81 

81 

81 

228 

228 

F s (kN) 

4 

8 

8 

8 

2 

4 

V w dg (nm.s -1 ) 

16 

16 

16 

16 

16 

16 

$ (%) 

1 

0 

0 

0 

1 

0 

Experiment n° 

13 

14 

15 

16 

17 


d 0 (/mi) 

228 

583 

583 

583 

583 


F s (kN) 

8 

8 

8 

8 

8 


Vwdg (nm.s -1 ) 

16 

1.6 

16 

160 

1600 


$ (%) 

0 

0 

0 

0 
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Experiment n° 

18 

19 

20 

21 

22 

23 

d 0 (/mi) 

21 

81 

42 

42 

42 

42 

F s (kN) 

0.1 

0.1 

0.1 

0.2 

0.6 

0.8 

Vwdg (nm.s -1 ) 

16 

16 

16 

16 

16 

16 

$ (%) 

11.65 

18.7 

15.1 

10.4 

^3 

«2 


TABLE I. Values of the experimental parameters associated 
with the different fracture surfaces analyzed here: bead diam¬ 
eter do prior to sintering, sintering force F s , wedge velocity 
Vwdg, porosity <f>. 


of crack propagation, y perpendicular to the mean crack 
plane, and z parallel to the mean crack front. 

The sample is then placed between the two jaws of 
a loading machine (home-made one or commercial elec¬ 
tromechanical machine). A home-made metallic wedge 
(semi-angle of 15°) is placed between the jaw and the 
specimen notch (Fig. Hf). It has been designed to min¬ 
imize friction dissipation by means of two rollers and to 
convert the compression along x into a tension along y 
[391140]. The sample is then compressed by lowering the 
upper jaw at a constant velocity. In the experiments re¬ 
ported here, the wedge speed V w d g varied from 1.6 nm.s -1 
to 1600 nm.s -1 . During each test, the force F(t) applied 
by the wedge to the specimen was monitored in real time 
by a static load cell. As soon as the wedge starts to push 
on the specimen, F increases. When F(t) reaches the 
critical loading (Fig. If)- the crack starts to propagate. 
Prior to fracture the load-displacement curve is linear, 
confirming the elastic behavior of the material. The crack 
velocity is observed to first increase (over ~ 20 mm), then 
to stabilize at a roughly constant value, ~ 100 times the 
wedge speed Wdg (over 40 to 60 mm), and to, finally, 
decrease over the last 20 mm. Our study focus on the 
interval where the velocity is constant and all measure¬ 
ments are done in that region. 


C. Post-mortem analysis of the crack surfaces 

Once the samples were broken into two halves, the 
morphology of the post-mortem crack surfaces were an- 




Wedge displacement (mm) 

b) 

FIG. 2. (Color online) a) Experimental setup sketch, b) Typ¬ 
ical force vs. displacement curve. 


alyzed. Firstly, the surfaces were visualized by means of 
a numerical microscope Keyence®. This microscope au¬ 
tomatically performs a vertical scan at different heights 
and reconstructs a well-focused image pixel by pixel over 
the whole topography range, irrespectively of the level of 
height variations. 

Topography profiles are also recorded by a Bruker me¬ 
chanical profilometer. Three profiles are scanned along 
the x-direction and 3 others along the ^-direction. Their 
lengths are respectively 55 mm along x, and 14 mm along 
z. Each profile has been positioned in the central part 
of the fracture surface away from the sample borders. 
They are located between x = 40 and 95 mm and at 
z = {—3, 0, 3} mm for the profiles in the x direction, and 
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between z — —7 and 7 mm and at x = {—40,60,80} 
mm for the profiles in the z direction. The out-of- 
plane and in-plane resolutions are Sy = 0.1 nm and 
Sx = 5z = 0.8 /im, which corresponds, per profile, to 
68750 points in the x direction and 16250 in the z direc¬ 
tion. 

For do = 583 /im, the fracture surfaces have been mea¬ 
sured using another profilometer due to the vertical dis¬ 
placement limitation of the Bruker profilometer [30] . Its 
out-of-plane and in-plane resolutions are Sy = 100 nm 
and Sx = Sz = 1 /im. To have a better statistic, ten 75 
mm long profiles have been measured in the x direction. 
Yet, in this case, due to the small number of beads in the 
thickness, no measurements were made in the z direction. 


III. RESULTS 

A general overview of the fracture surfaces obtained by 
changing the bead size, the porosity and the crack veloc¬ 
ity is first given. In most of the cases, the propagation 
is inter-granular. We focus on this case, investigating 
the effects of the microstructure size and porosity on the 
fracture surfaces. 


A. Fractography 

Figure [3] presents the microscope visualizations of typi¬ 
cal fracture surfaces, for different bead diameters, porosi¬ 
ties and wedge velocities. The crack propagates from left 
to right. 

The snapshots a,b,c correspond to a null porosity and 
to samples broken at the same wedge velocity, but with 
different bead diameters. The particles surfaces are ob¬ 
served to be completely deformed resulting in sharp edges 
and corners, in agreement with the zero value of the 
porosity. Also, the surfaces are found to be made of 
smooth facets. Such a faceted morphology is characteris¬ 
tic of a brittle inter-granular fracture, with a crack grow¬ 
ing along the grain-grain interfaces. Note that the facets 
are slightly more elongated along x than along z (visi¬ 
ble on Fig. IF This anisotropy is due to the sintering 
process. In fact, the compression was applied along z, 
or in a granular packing, the two components a xx and 
a zz of the stress tensor are classically related m via 
a xx = Kcr zz where the Janssen constant K is smaller 
than 1, typically around 0.6-0.8. During sintering, the 
beads, hence, contract more along z than along x. The 
snapshot (d) of Fig. [3j corresponding to a smaller wedge 
speed F^dg = 1.6 nm/s presents a qualitatively differ¬ 
ent morphology: The facets are observed to be blurred 
with a multitude of small-scale fragments. They betray a 
intra-granular propagation mode, with a crack propagat¬ 
ing throughout the grains, without necessary following 
the grain-grain interfaces. 

Figure [3^, f shows microscope visualizations of fracture 
surfaces of sintered samples realized with lower sintering 


force F s . At low F s (Fig. the beads globally keep 
a spheric shape. The contact areas between adjacent 
beads is small and are disc-shaped. As a result, an im¬ 
portant porosity is found in such samples. When F s gets 
larger, the contact surfaces grow and intersect, result¬ 
ing in the apparition of sharp edges of increasing length 
(Fig. [3]f), while the undeformed spherical parts of the 
particles shrink, so that the volume of the pores, delim¬ 
ited by these parts, decreases. When F s is large enough 
to close all the pores (Fig. [3^) , the undeformed parts of 
the particles disappear and the corners are sharp. 

Figure [4] presents typical topographical profiles H(z) 
as measured along the z direction on the fracture surfaces 
observed in Fig. [3j In the Fig. [4}i,b,c associated with 
inter-granular fracture, the corrugations in the y and z di¬ 
rections are of the same order, that of the bead size. Con¬ 
versely, the roughness observed for intra-granular mode 
(Fig. §1) appears much flatter, and exhibits plateaus 
much longer than the bead diameter. These plateaus 
correspond to zones where the crack has cut through¬ 
out the grains, without being perturbed by the inter¬ 
faces. Finally, Fig. [4^,f show typical profiles at finite 
porosity. Note the increase in roughness with amplitudes 
larger than the bead diameter do. The shift from inter- 
to intra-granular fractures occurs for fracture velocity of 
1.6 nm/s. In the present study, we focus on the inter¬ 
granular case, hence on experiments corresponding to 
Vwdg >1-6 nm/s. Works on the transition between intra- 
and inter-granular propagation and on the intra-granular 
fracturation are in progress. 


B. Effect of the microstructure length-scale 

We first restrict the analysis to the effect of the mi¬ 
crostructure length-scale and, hence, only consider the 
specimens (1-17 in Tab. [T|) with very low porosity (< 
0.5%). To quantify the spatial distribution of crack 
roughness, we computed the structure functions [42] 
S z (Az) along z, and S x (Ax) along x, defined by: 

S z (Az) = ((H(z + Az)-H(z)) 2 ), (1) 

S x (Ax) = ((H(x + Ax)-H(x)) 2 ), (2) 

where the operator () refers to an average over all possible 
positions 2 and x, respectively. 

Figure [5] displays in a log-log plot the structure func¬ 
tions S x measured along x for different bead diameters. 
Two scaling regimes can be distinguished. At large scale 
(Ax 10/im), S x slowly increases with Ax, with a pref¬ 
actor increasing with do. This betrays the fact that the 
roughness scales with do. At small scales (Ax <C 10/im), 
all curves collapse, the roughness is in this regime weakly 
dependent on do. 

The diameter do appears to be the only length-scale 
in the problem. Therefore, we made dimensionless the 
structure functions by making Ax —> Ax/do and S x —>> 
S x /do (Inset of Fig. pi). A good collapse is obtained, con- 
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a) Experiment n°2, do=21/im, b) Experiment n°13, do=228/im, c) Experiment n°16, do=583/im, 

V w dg=16 nm.s -1 , F s = 8 t, 4> = 0% V w dg=16 nm.s -1 , F s = 8 t, 4> = 0% V w d g =16 nm.s -1 , E s = 8 t, 4> = 0% 



d) Experiment n°14, do=583/im, e) Experiment n°21, do=42/im, f) Experiment n°3, do=42/im, 
V w dg=l-6 nm.s -1 , E s = 8 t, <t> = 0% V w d g =16 nm.s -1 , F s =0.1 t, <t> = 15% V w d g =16 nm.s -1 , F s =l 1. <f> — 1.45% 


FIG. 3. Microscope visualization of wedge-splitted fractured surfaces of the sintered-PS-beads samples. The bead diameter do, 
the wedge speed E w d g , the sintering load F s and the corresponding porosity 4> are indicated beneath each snapshot. 


firming that do is the length-scale governing the spatial 
distribution of the roughness. 

Nevertheless, the collapse observed in the inset of Fig. 
[5] is not perfect and a clear scattering remains visible. We 
applied the global minimization technique described in 
[43] to reduce the dispersion of the data. For each curve, 
an optimal length £ is obtained. The values can be found 
in Table [TT] for each sample. Fig. [ 6 ^ 1 , b shows the structure 
functions in the x and z directions made normalized by 
£ x and i z . As can be seen, the sole parameters £ x and i z 
are sufficient to achieve a very good collapse. We checked 
that these lengths correspond to the microstructure sizes 
in both directions (as e.g. observed in fig. [3|. In other 
words, the slight differences observed between £ x , i z and 
do result from the sintering process. 

Moreover, two power-law regimes, characterized by two 
different exponents, are observed. To determine these ex¬ 
ponents, we performed linear fits of the data displayed in 
Fig. § The roughness exponents £+’ and £+’ , for 
length scales above and below do, are defined by the re- 


Experiment n° 

1 

2 

3 

4 

5 

6 

d 0 (/xn 1 ) 

21 

21 

42 

42 

42 

81 

$ (%) 

1 

0 

1 

0 

0 

1 

4 (/mi) 

19.8 

19.6 

19.6 

39.9 

44 

79.1 

4 (/mi) 

20.2 

21.8 

26.4 

44.8 

44.9 

81.7 

Experiment n° 

7 

8 

9 

10 

11 

12 

d 0 (/im)l 

81 

81 

81 

81 

228 

228 

$ (%) 

1 

0 

0 

0 

1 

0 

4 (/mi) 

70.2 

80.3 

68.4 

76 

203.3 

187.4 

4 (/mi) 

74.8 

82.9 

81.7 

83.7 

238.7 

229.2 

Experiment n° 

13 

14 

15 

16 

17 


d 0 (/mi) 

228 

X 

583 

583 

583 


$ (%) 

0 

X 

0 

0 

0 


4 (/mi) 

166 

X 

X 

X 

X 


4 (/in 1 ) 

192.9 

X 

381 

486.5 

440.8 



TABLE II. Optimal lengths £ x and £ z along the x and z di¬ 
rections used to collapse (in Fig. [6j the structure functions. 
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_1 Typical 

a) Experiment n°2, d Q =21 M m, V wdq =16 nm.s , 0=0% ^ad size 




b) Experiment n°13, d n =228p m, V wdg = 16 nm.s , 0=0% 




c) Experiment n°16, d =583pm,V =16 nm.s \ 0=1 



d) Experiment n°14, d n =583p m, v wdq = 1-6 nm.s , 0=0' 



^ Op 


e) Experiment n°21, d 0 =42p m, V wdg = 16 nm.s ,0=15% 




f) Experiment n°3, d Q =42 jlx m, V wdg =16 nm.s , 0= 




FIG. 4. Experimental height profiles along the crack propa¬ 
gation direction of the same samples than in Fig[3] On the 
right, the bead shape is given as reference, using the same 
abscissa and ordinate scales as in the plots. Since these scales 
are different, the bead appears distorted. 




a) 



10 

Az/t z 

b) 


FIG. 6. (Color online) Structures functions normalized by the 
optimal lengths i x and i z . a) and b) are respectively for the 
measurements performed along x (propagation direction) and 
z (front crack). Both axes are logarithmic in the main panels, 
and semi-logarithmic in inset. The straight lines in the main 
panels are power-law fits with exponents reported in Tab. |III| 


Direction 

X 

z 

Small scales : 

0.86T0.04 

0.81T0.05 

Large scales : C + 

0.27T0.02 

0.35T0.01 


TABLE III. Roughness exponents obtained by linear fits per¬ 
formed on the curves of Fig. [6] 


FIG. 5. Structure functions S x of the inter-granular exper¬ 
iments, along the propagation direction at 4> = 0%. Each 
symbol corresponds to a given value of do. Inner window : 
Same distributions made dimensionless by de¬ 


lations S x oc Ax 2< * x and S z 


oc A z 2( * z . 


Table III 


gives 


the roughness exponent values along the propagation and 
front crack directions. In the small scales regime, the ex¬ 
ponents are found close to 0.8 in both directions. At 
larger scales, the exponents £ + are functions of the di¬ 
rections of measurements, the roughness along the crack 
being more important than the one in the propagation 
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direction. Although being relatively small, the insets of 
Fig. © clearly prove that they are finite and exclude a 
logarithmic behavior. 

C. Effect of the porosity 

We now investigate the influence of the porosity <F. 
Fig. [7] shows the dimensionless structure functions. Con¬ 
trary to the previous case, the data do not collapse. First, 
the amplitude of the structure functions is found to in¬ 
crease with <f> at all scales. This is compatible with the 
observations of rougher profiles as the porosity increases 
(Fig. §. At small length-scales (Ax/do < 1 ), the expo¬ 
nents are similar for all porosities. A linear regression of 
the data gives a roughness exponent = 0.88 db 0.04, 
very close to the one otained for <£> = 0. Yet, for 
Ax/do > 1 , the exponents rise with the porosity. 

In Fig. [ 8 j we analyze more quantitatively the effect of 
<f> on the roughness. The roughness amplitudes, defined 
as the value of the structure functions at Ax/do = 1, are 
observed to increase linearly with <F (Fig. [ 8 ^ 1 ). Fig. [ 8)3 
shows the evolution of the large scale roughness expo¬ 
nents C+ and C+ as a function of <F. As expected, when 
<f> tends toward 0 , the values obtained in the no porosity 
case (Table. are recovered. They increase with <F, up 
to values (+ ~ £+ ~ 0.48 for <F = 20%. No significant 
anisotropy is observed on the roughness amplitude (Fig. 
§L). Conversely the scaling anisotropy, i.e. a difference 
between £+ and £+, observed for <F = 0 is recovered. 
However, it decreases as <F increases. 

To quantify the scaling anisotropy evolution, we plot 
the ratio £+ /£+ as a function of <£> in Fig. [ 9 ] It decreases 
from C+/C+ = 1.3 ± 0.05 to 1 as $ goes from 0 to 10%. 
Above, isotropic scaling is obtained. 

IV. DISCUSSION 

A. Comparison with other experimental results 

The roughness value £“ measured at small scales is 
observed to be independent of the bead size, the direc¬ 
tion measurement, and of the porosity. It is compatible 
with the C = 0.8 value reported in the literature for a 
variety materials Here, £ is observed at a scale 

smaller than the bead size. This small scale regime is not 
strictly speaking a self-affine regime: The value can 
be interpreted as the measured roughness exponent of a 
piecewise linear profile analyzed at scales smaller than 
that of its segments (see m for related analysis). 

Regarding the roughness properties at large scales, 
it is interesting to compare our results with what has 
been reported for ceramics obtained by sintering oxide 
glass beads [28] , The roughness amplitude (defined as 
yjS(Sr = do)) is found to increase linearly with <F as ob¬ 
served in ||28] for glass ceramics. Conversely, our data 
converge toward a finite amplitude as <F 0 , while, in 



Ax/d 0 


FIG. 7. Dimensionless structure functions along the propaga¬ 
tion direction for different porosities <£>. 


glass ceramics, the amplitude was vanishing. The main 
difference is that, in Ref. [28] 5 decreasing <f> induces a 
transition from an inter to intra granular propagation 
mode, this is not the case here. 

The values of the roughness exponents are also, within 
the errorbars, similar to the values ( ~ 0.4 —0.5 reported 
in these glass ceramics [28] and in sandstone [27] l30j. Our 
material is hence more representative of conventional het¬ 
erogeneous brittle materials such as ceramics and some 
rocks. As a bonus, the smaller errorbars and the larger 
range of <F explored have evidenced the increase of ( with 
<F. They have also revealed a scaling anisotropy in ab¬ 
sence of any porosity, which decreases as <F increases and 
disappears for <f> > 10 %. 


B. Comparison with theoretical results 

Since the early 1990s, a large number of theoretical 
work has focus on the spatial distribution of fracture sur¬ 
face roughness in brittle heterogeneous materials. They 
can be classified into two categories: (i) elastic string 
models that consider the crack front as an elastic line 
propagating through randomly distributed obstacles la- 
mmm-, and (ii) random lattice models that model the 
material by a network of fuses, springs or beams with ran¬ 
domly distributed breakdowns thresholds Eng Dana. 
The main difference between the predictions of these two 
classes of models is that models of type (i) naturally lead 
to anisotropic surfaces, where the direction of front prop¬ 
agation plays a specific role, while model of type (ii) lead 
to isotropic surfaces m • 


























a) 



b) 


FIG. 8. a) Normalized roughness amplitudes in the x and 
z directions \/S(do)/do versus porosity <f>, at the position 
Ax/do = 1. b) Large scale roughness exponents vs. <f>. 
Markers match the experiments of Fig. [7] The crosses corre¬ 
spond to the averaged exponent values of Table |III| 


By essence, elastic line models address situations of 
nominally brittle fracture, with nonporous materials. 
Hence, they are the relevant theoretical framework to 
discuss the observations reported in Sec. |IIIB| In partic¬ 
ular, elastic string approaches lead to: 


a) The length scale of the microstructural disorder to 
be the single relevant length-scale for the structure 
function; 


b) Anisotropic surfaces, where the direction of front 
propagation plays a specific role. 

Both these predictions are in agreement with the exper¬ 
imental observations (Figs. [6]and|8|. 



FIG. 9. Scaling anisotropy defined as (£/(+ as a function of 
the porosity <f>. 


Conversely, the values of the roughness exponents 
C+ « 0.27 and £+ ^0.35 observed here are clearly dis¬ 
tinct from the values predicted theoretically. In particu¬ 
lar, the most refined models mm attempted to derive 
rigorously the equation of motion of the elastic string 
from Linear Elastic Fracture Mechanics. They yield ei¬ 
ther logarithmic scaling 12 or ( x = 0.50 ± 0.05 and 
( z ~ 0.385 ±0.05 [14] according to the disorder intro¬ 
duced in the equation. The values reported here are 
clearly different, which indicates that some key ingre¬ 
dients are missing in the models. In particular, no model 
takes into account the T-stress influence. The T-stress 
here has been determined m both by finite elements 
simulations and digital image correlation [48] . As shown 
in Fig. [lO] the T-stress is positive, so that the crack 
path should be unstable toward any perturbation of its 
rectilinear path [49]. In other words, additional corre¬ 
lation are anticipated to result from a positive T-stress, 
and hence a larger value for the roughness exponent mea¬ 
sured along x. If this scenario is correct, it may partly 
explain the departure from the logarithmic behavior. 

As mentioned above, in elastic line models, the dimen¬ 
sionless curves S x /d^w s Ax/do and S z /do vs Az/do are 
anticipated to be universal. But, we observed a scaling 
function of <£> (See Fig. [7]). This indicates that elastic 
line approaches stop being relevant as the porosity be¬ 
comes finite. This is thought to be due to the presence of 
microcracks forming at the pores, and making the frac¬ 
ture propagation mode shifting gradually from nominally 
to quasi-brittle. The anisotropy decrease with <£> is con¬ 
sistent with this scenario: As microcracking develops, 
the fracture surface stop being the trace left by a single 
propagating line in an otherwise intact material, but re¬ 
sult from the coalescence of multiple microcracks instead. 
Random lattice models include such processes and, may 
be the relevant framework to describe this quasi-brittle 
regime. In this context, it is interesting to note that 
the roughness exponents « 0.45 measured here 
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FIG. 10. (Color online) T-stress as a function of the crack tip 
position xtip obtained by classical finite element simulations 
and digital image correlation. 


on isotropic fracture surfaces, i.e. above = 10%, are 
consistent with the values Q oc = 0.42, (i oc = 0.5 and 
Cioc = 0-48 observed respectively in 3D random fuse m, 
spring m or beam m networks. 

In the above scenario, porosity provides a tunable pa¬ 
rameter to go gradually from nominally brittle fracture, 
where the disorder effect is to distort the front propaga¬ 
tion, to quasi-brittle fracture where the fracture surfaces 
emerge from a percolating path throughout the micro¬ 
crack cloud. At present, the theoretical descriptions of 
these two situations belongs to two distinct realms : Elas¬ 
tic line models for the former and random lattice models 
for the latter. Unifying these two frameworks represent 
an important challenge for future investigations (see [36] 
for a recent thrust in this direction). 

V. CONCLUSION 

In this paper, we have used home-made cemented 
grains materials to investigate the influence of the mi¬ 
crostructure size and porosity on the fracture surfaces. 
The following observations have been made: 

• The roughness displays scale invariant morpholog¬ 


ical features that depends on both the grain size 
and porosity. 

• At zero porosity, the structure functions measured 
along a direction and normalized by the grain size 
collapse onto a single master curve. 

• Scaling anisotropy is well pronounced for <£> = 0, 
but decreases with and vanishes for 4> > 10%. 

• The values of the roughness exponents increase 

with 4>, from £+ ~ 0.27 and £+ ~ 0.35 at 4> = 0% 
to C ~ 0-48 at 4 > = 20 %. 

These observations are compatible with previous stud¬ 
ies considering other sintered glass beads materials [28] 
or sandstone Hang. But in these previous studies, the 
precision was not sufficient to detect any clear varia¬ 
tion with <£>. However, the values of the roughness ex¬ 
ponents C+ « 0.27 and £+ « 0.35 observed, at low 
porosity, are clearly distinct from the values predicted 
theoretically : additional investigations are necessary to 
take into account some missing ingredients like the T- 
stress. Conversely, for <£> > 10%, the roughness expo¬ 
nents C+ = Ct ~ 0.45 are consistent with the values 
(i oc between 0.4 and 0.5, obtained with random networks 
mils iso]. Modulating 4> has permitted, thus, to mod¬ 
ify quantitatively the value of the roughness exponents. 
Note that, in conventional interface growth problems, the 
roughness exponent value is characteristics of a univer¬ 
sality class. Their continuous evolution with 4> observed 
herein may pose a rather severe test for current and fu¬ 
ture competing models of heterogeneous fracture. The 
T-stress component acting on the fracture is controlled 
by the macroscopic shape of the porous samples; future 
works will consider the effect of the T-stress on the rough¬ 
ness. 
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